Theory of a large thermoeffect in superconductors doped with magnetic impurities 
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We argue that parametricaUy strong enhancement of a thermoelectric current can be observed in 
conventional superconductors doped by magnetic impurities. This effect is caused by violation of 
the symmetry between electron-like and hole-like excitations due to formation of subgap Andreev 
bound states in the vicinity of magnetic impurities. We develop a quantitative theory of this effect 
and demonstrate that it can be detected in modern experiments. 
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Application of an electric field £^ to a normal conductor 
with Drude conductivity (Jat yields an electric current 
j = (tnE across this conductor. A similar effect can be 
produced by a temperature gradient VT. In this case the 
current j induced in a sample takes the form j = a^VT^ 
where a^ ~ {cfn /e){T /ep) is thermoelectric coefficient 
and ep is the Fermi energy. The latter simple equations 
illustrate the essence of the so-called thermoelectric effect 
in normal metals. 

If a metal becomes superconducting, the situation 
changes significantly. On one hand, the electric field can- 
not anymore penetrate into a superconductor and, hence, 
the Drude contribution to the current is absent in this 
case. On the other hand, a supercurrent js can now be in- 
duced in the sample without any electric field. It follows 
immediately that by applying a temperature gradient to 
a uniform superconductor one would not be able to in- 
duce any current since thermal current would be exactly 
compensated by the supercurrent js = —aVT^ where 
a defines thermoelectric coefficient in a superconducting 
state. Ginzburg [l|, |2| demonstrated that no such com- 
pensation generally occurs in non-uniform superconduc- 
tors which opens a possibility to experimentally detect 
thermoelectric current in such structures. Several exper- 
iments with bimetallic superconducting rings (see Fig. 
[T]) have been performed [3-55 which indeed revealed the 
presence of thermoelectric magnetic flux in such rings. 
However, both the magnitude of the effect and its tem- 
perature dependence turned out be in a strong disagree- 
ment with available theoretical predictions [6[. Quite 
surprisingly, the magnitude of the thermoeffect detected 
in these experiments exceeded theoretical estimates by 
several orders of magnitude. Subsequently, a good agree- 
ment between theory and experiment j7] was claimed, but 
this report remained largely unnoticed. In any case, no 
convincing explanation of the discrepancy between exper- 
iments [3|-l5| and theory [6| was offered and the paradox 
remains unresolved until now [S^. 

In this Letter we are not aiming at directly resolving 
this long standing paradox. Rather our primary goal is to 
identify the conditions under which thermoelectric cur- 



FIG. 1: A ring formed by two different superconductors with 
contacts maintained at different temperatures Ta and T^. 



rents in superconductors can be significantly enhanced. 
In the normal state contributions to the thermoelectric 
coefficient a^ from electron-like and hole-like excitations 
are of the opposite sign and almost cancel each other. A 
similar situation occurs in conventional superconductors 
where the thermoelectric coefficient a also remains small 
[6] and monotonously decreases with T below the critical 
temperature Tc. On the other hand, in unconventional 
superconductors impurity scattering may lead to much 
larger values of a due to formation of quasi-bound An- 
dreev states near impurities which yield high asymmetry 
between electron and hole scattering rates [9, 10]. 

Here we will demonstrate that "giant" thermoeffect 
can also be expected in conventional superconductors 
doped by magnetic impurities. Also in this case Andreev 
bound states are formed near such impurities |llMl3J 
thereby explicitly breaking the symmetry between elec- 
tron and holes [M, JJ.] • We argue that this feature may 
cause parametricaUy strong enhancement of the thermo- 
effect in such systems 



a/aN{Tc) ^ Pf^ > 1, 



(1) 



where pf = tuvf is the Fermi momentum and i is the 
electron elastic mean free path in the absence of mag- 
netic impurities. This formula remains valid in the most 
relevant diffusive limit £ ^ vf/Tc and at an "optimal" 
concentration of magnetic impurities nimp roughly equal 



to one half of the critical one (see below). Eq. ([T]) pre- 
dicts possible enhancement of the thermoeffect in super- 
conductors with magnetic impurities by several orders of 
magnitude as compared to that in the normal state at 

Quasidassical formalism and impurity self- averaging. 
In what follows we will consider a superconductor which 
contains both non-magnetic and magnetic impurities. 
Our analysis is based on the quasidassical formalism of 
nonequilibrium Green-Keldysh matrix functions g obey- 
ing the Eilenberger equations [16] 

-ivFVg{pF,r,£,t) = [(l-t,g{pF,r,e,t)\, f = 1. (2) 

The check symbol denotes 4x4 Keldysh matrices 

^=[o xa)^ X = g,n,^, (3) 

with blocks X^'^'^ being 2x2 matrices in the Nambu 
space. The matrix Cl has the standard structure 



n'' = n' 



e A 

-A -£ 



n 



K 



0, 



(4) 



where e is the quasiparticle energy, A is the BCS order 
parameter which is chosen real further below. 

Scattering of electrons on impurities is accounted for 
by the self-energy matrix £ which can be expressed in 
the form 



-iV^ 



^ vf/{2£). 



(5) 



Here the first term describes the effect of non-magnetic 
isotropic impurities while the second term £m is responsi- 
ble for electron scattering on randomly distributed mag- 
netic impurities ^13j 



Zjm — 



27TNi 



rM 



ui + T3U2] ^ + i (<?)) + 

+ {[ui-f3U2]-' + i{g)y^}, (6) 



where A^o is the electron density of states per spin direc- 
tion at the Fermi level, Ui^2 are dimensionless parameters 
characterizing the impurity scattering potential and fa is 
Pauli matrix in the Nambu space. Averaging over the 
Fermi surface is denoted by angular brackets (•••). Note 
that within the Born approximation the self-energy (J6j) 
just reduces to the well known Abrikosov-Gor'kov result 
[17]. Unfortunately this approximation is insufficient for 
our present purposes since it does not allow to account 
for impurity Andreev bound states (impurity bands) and 
the electron-hole asymmetry. For this reason in what fol- 
lows we will go beyond Born approximation and employ 
a more general expression for the self-energy (|6]). 

Finally, the current density j is defined with the aid of 
the standard relation 



j{r,t) = ^ / de{vFSp[T3g^{pF,r,s,t)]) . 



(7) 



Electron-hole asymmetry and the density of states. It is 
well known that two subgap Andreev bound states with 
energies 



SB = ±/3A, P^ = 



(1 



,2\2 



(1 



uir+4ur 



(8) 



are localized near each magnetic impurity in a super- 
conductor [ll|, ll2|. Similarly to the case of uncon- 



ventional superconductors with non-magnetic impurities 
[9|, llO| these Andreev bound states yield different scat- 
tering rates for electrons and holes and, hence, break the 
electron-hole symmetry in our system thereby causing 
strong enhancement of the thermoeffect. 

Consider the retarded part of the self-energy I] ([5]). It 
can be written in the form 



^i? 



So^ 



yR 

^F+ 



^f 



yR 



yR 



(9) 



where non- vanishing diagonal part T,q exphcitly accounts 
for asymmetry between electrons and holes [15]. Substi- 
tuting the retarded Green function matrix 



A2 



e A 
-A -s 



(10) 



into Eqs. (J5]), (|6]) we evaluate T^q as well as the energy 
resolved superconducting density of states u{e) normal- 
ized to its normal state value. Introducing the parameter 

e = sA/A we get 



s^(£) 



0^ 



v{e) = Re 



A2 = y; 



/32A2' 

A2 + ir + iri- 



Ve 



?2 



2 - A2 
A2 



/32A2' 



(11) 

(12) 



where the parameter e is fixed by the relation [ll|, [12 1 

(13) 



iV 



V^ 



eye 
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The scattering parameters ro,i,2 have the dimension of 
rates being proportional to the concentration of magnetic 
impurities nimp- They read 



To 



Ti 



wi(l 



ul) 



ttNo (l + uj- u^y + 4^2 ' 



ap(l 



ul){ul - ul) 



ttNo (1 + u1 - ul)^ + 4ul 



T. = 2- 



nNo (1 + wf - w^)2 + 4ui 



(14) 
(15) 
(16) 



Note that the parameters £, J^q and u{e) remain insensi- 
tive to the electron scattering rate on non-magnetic im- 
purities r since such scattering does not produce any 



Uiu^p = O.IttNoTcO ;■ 

riiuip = O-^ttNqTco ; 
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FIG. 2: (Color online) Energy resolved density of states 
iy{e) = ^{—£) in a superconductor doped by magnetic im- 
purities. Tco is the critical temperature of an undoped super- 
conductor. 



pair-breaking effect in bulk conventional superconduc- 
tors. On the contrary, scattering on magnetic impurities 
may strongly modify these parameters. For illustration, 
the density of states ^{e) is depicted in Fig. [2] at £ > 
and different values nimp- With increasing nimp Andreev 
levels get broadened forming two impurity bands respec- 
tively at positive and negative energies. Further increase 
of nimp yields even broader bands which eventually merge 
with continuum (overgap) states. 

Thermoeffect enhancement by magnetic impurities. 
We are now prepared to evaluate the thermoelectric co- 
efficient a. In doing so we will essentially follow the qua- 
siclassical linear response theory initially formulated in 
Ref. 118| for the analysis of thermal conductivity in un- 
conventional superconductors. This approach allows to 
recover the dominating contribution to the thermoelec- 
tric coefficient a which originates from the electron-hole 
asymmetry. Employing Eqs. (|2j) and proceeding along 
the lines with Ref. 



181 we evaluate the correction to the 
Keldysh Green function Sg^ ex VpVT [19|. Combining 
the resulting expression with Eq. ([7]) we obtain 



a ■■ 



He) 



eNpvj, r^ J^{e)de 
' 12T2 J_^cosh\e/2Ty 



Re./A2 



[lm^^{e)y 



(17) 
(18) 



Eqs. ([I71)-([l8l) - together with Eqs. ([II])-([l6D - con- 
stitute the central result of this work which accounts 
for "giant" thermoeffect in superconductors doped by 
magnetic impurities. In the most relevant case of dif- 
fusive superconductors with F ^ Tc Eq. (fTSj) reduces to 
J='{s) = u{s)lmJ:^{e)/T'^, i.e. a ex l/F^ in this limit. 

At small magnetic impurity concentrations r2 <C A 
the impurity band is restricted to subgap energies e < A, 
cf. Fig. [21 At energies within the impurity band one has 



He) 



To 



2r2AVr^^ 



/3A)2 



4r2(i-/?2)(Av/r^^ + r)2 



(19) 



Substituting this expression into Eq. (fT7|) . integrating 
over all impurity band energies and taking the limit F2 <C 



nimp = O.bTTNoTcO 

nimp = I.OttNqTco ■ 

nimp = 2.^7tNoTcO ■ 
nimp = S.OTTiVo^cO 
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FIG. 3: (Color online) Thermoelectric coefficient as a function 
of temperature and magnetic impurity concentration. Scat- 
tering parameters ui = U2 = 0.5 and the scattering rate 
r = lOTco are the same for both panels. 



A,T^/A, we arrive at the subgap contribution to a: 



^sg 



eNpvj, cosh~^(/3A/2T)Fov/2r^A^/^ 
9T' (l-/32)i/4(Ayr^^ + F)2' 



(20) 



i.e. asg ex n-^ at small concentrations of magnetic im- 
purities. Assuming that the impurity band is located at 

e ~ A/2 {P ~ 0.5) and setting T ~ A ~ To, 1,2 we get 



-eNovj.T/r'^ - aNPF^' 



(21) 



This estimate demonstrates that asg may strongly exceed 
the thermoelectric coefficient in the normal state. 

Additional contribution to a is provided by overgap 
energies. For small values nimp we can use the standard 
BCS expression for the density of states and derive 

2ror2(l - /3^)e^A^[e^ - /3^A^]-^ 



He) 



[T2e^ + r(£2 - /32A2) + Fi(£:2 - A2)] 



2 • 



(22) 



Combining Eqs. (|22]) and (p!7|) . in a realistic limit F ^ 
ri^2 and for T '^ A we recover the contribution to a from 
overgap energies 



Tor 



OLo 



QJ- 2 



-eA^o4^p 



ror2 . 

OiN—r^-PF^' 



(23) 



In the optimal case T ~ A ~ ^0,1,2, we find ao 



where the latter quantity obeys Eq. (|2T]) . Hence, also for 
a = asg + aog we recover the estimate ([T]). 

At temperatures close to Tc the value a can be evalu- 
ated analytically at any concentration of impurities. In 
this limit one can set u{e) = 1 and obtain 



Tie) 



2ror2(i - /?^)£^A^ 
(e2 + rl)2(r + ri + r2)2 



(24) 



which yields 



eiVo4 ro(l - /3^)A^ 
' 67rT2 (r + Ti 



■r2)2 



s 



2itT 



(25) 



where S{x) = [xijj' {x + 1/2)] and '^(x) is the digamma 
function. 

The results of numerical evaluation of a as a func- 
tion of both temperature and impurity concentration are 
displayed in Fig. [3l We observe that the thermoelec- 
tric coefficient of a diffusive superconductor achieves its 
maximum value at temperatures T ~ Tc/2 and nimp ap- 
proximately equal to one-half of the critical concentration 
at which superconductivity gets fully suppressed. This 
maximum value can be estimated as 



max \a\ 

J- i'^T'iniD 



0.05 



eNovjrTco 



0.2eiVoTeor 



(26) 



Combining the expression for a^ ^ {(Jn /^)iT /ep) with 
Eq. (|26|) we arrive at the estimate ([!]) which demon- 
strates that enhancement of the thermoeffect is stronger 
in cleaner superconductors. At the borderline of applica- 
bility of Eq. ([T]) £ ~ vp/Tc we obtain \a\ r^ a^ je^ which 
appears to define the absolute maximum value of a in 
conventional superconductors doped by magnetic impu- 
rities. 

It is interesting to point out that the presence of 
electron-hole asymmetry in such superconductors was 
also predicted to yield anomalously large photovoltaic 
effect [15]. Despite clear similarity between the mod- 
els the effect [15] is substantially different from one ana- 
lyzed here. Indeed, while no voltage occurs in the system 
within the linear response to a temperature gradient [19], 
a non-zero nonequilibrium voltage is induced as a second 
order response to an external electromagnetic field [15]. 
Hence, thermal heating of the system considered here is 
physically not equivalent to that produced by an external 
ac field. 

Bimetallic superconducting rings and TEB. Finally let 
us briefly discuss the possibility to experimentally detect 
"giant" thermoeffect predicted here. One way to do so 
would be to perform an experiment with bimetallic su- 
perconducting rings [3-5] as shown in Fig. [H Provided 
superconducting contacts are kept at different tempera- 
tures Ta and T5, thermoelectric current will be induced 
inside the ring and the corresponding magnetic flux ^ 
can be measured. The magnitude of this flux reads 



$0 



4e 



n 



— = -. [A?(T)ai(T) - \l{T)a2{T)]dT, (27) 



Ta 



where <l>o = 7rc/e is flux quantum, ai^2 and Ai^2 define 
respectively thermoelectric coefficients and the values of 
London penetration depth for two superconductors. For 
simplicity we may assume ai ^ 0^2 and neglect the sec- 
ond term in Eq. (|27|) . Employing Eq. ([26]) together 




nimp = O.SttTVoTcg ■ 

nimp = I.OttTVoTco ■ 

nimp = 2.57r7VoTco ■ 

nimp = S.OttTVoTcG 



FIG. 4: (Color online) Temperature dependence of the term 
— 4eaA^Tco/c^. Different curves correspond to different values 
nimp- The parameters 1^1,2 and F are the same as in Fig. O 



with the standard expression for the London penetration 
depth in diffusive superconductors at T = we arrive at 
a conservative estimate for the thermally induced flux 



M.0.01^ 

$0 r 



Ta\ 



T>T, 



cO- 



(28) 



In Fig. |4]we display the temperature dependence of 
the combination A^(T)a(T) at different concentrations 
of magnetic impurities. Induced thermoflux <I> (normal- 
ized to <I>o) equals to the area under the corresponding 
curve between Ta and T5. For reasonably clean supercon- 
ductors typical values of ^ may easily reach $ ^ 10~^<I>o. 

Another way to experimentally test our predictions 
would be to employ a novel type of zero-biased thermo- 
electric bolometer (TEB) [20|. This TEB consists of a 
superconducting absorber attached to normal and su- 
perconducting electrodes via tunnel junctions (SIN and 
SIS' junctions). Incoming photons excite quasiparticles 
in the absorber. Strong charge imbalance between ex- 
cited quasiparticles and quasiholes can be expected pro- 
vided this absorber is formed by a superconductor doped 
by magnetic impurities. Temperature gradient across the 
superconductor will occur due to permanent escape of ex- 
cited quasiparticles from the "cold" end of the absorber 
through SIN junction, whereas no such escape would be 
possible in its "hot" end attached to SIS' junction. As 
a result, permanent thermoelectric current will flow in 
the absorber creating "giant" thermoelectric voltage re- 
sponse which can be detected experimentally. 

In summary, we have demonstrated that "giant" ther- 
moeffect might occur in conventional superconductors 
doped by magnetic impurities. This effect is well in the 
measurable range and can be detected in modern exper- 
iments. 
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